Softer than normal, but not as soft as one might think: Spontaneous flux lattices in 

ferromagnetic spin-triplet superconductors 
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A theory is developed for the spontaneous vortex lattice that is expected to occur in the fer- 
romagnetic superconductors ZrZn2, UGe2, and URhGe, where the superconductivity is likely of 
spin-triplet nature. The long-wavelength fluctuations of this spontaneous flux lattice are predicted 
to be huge compared to those of a conventional flux lattice, and to be the same as those for spin- 
singlet ferromagnetic superconductors. It is shown that these fluctuations lead to unambiguous 
experimental signatures which may provide the easiest way to observe the spontaneous flux lattice. 

PACS numbers: 



Conventional wisdom holds that ferromagnetism and 
superconductivity cannot coexist 111 |2|. However, re- 
cent experimental studies of UGe2 [21 L| , URhGe [3, and 
ZrZn2 0| , backed by band structure calculations 0, @ , 
have demonstrated coexistence of both types of order in 
the same electron band in these materials. Theoretical 
considerations suggest that the superconducting order 
parameter is of the spin-triplet non-unitary type [|| ■ 

A state which displays coexistence of ferromagnetism 
and superconductivity is expected to have many unusual 
features, among them a spontaneous vortex or flux lat- 
tice. The internal magnetic field generated by the spon- 
taneous magnetization makes topological excitations, viz. 
vortices, in the superconducting order parameter energet- 
ically favorable |lfj. Unlike the well-known Abrikosov 
flux lattice 

.11]/ 

this spontaneous flux lattice state re- 
quires no external magnetic field. For spin-singlet su- 
perconductors, such spontaneous flux lattices have been 
proposed and theoretically studied previously 0,E3>E|- 

This Letter addresses the heretofore unstudied prob- 
lem of spontaneous flux lattices in spin-triplet p-wave 
superconductors. One might expect this problem to be 
much more complicated than the spin-singlet case, since 
the order parameter is much more complicated, which 
should lead to many more soft modes. One of our chief 
conclusions is that, surprisingly, this is not the case. 
Rather, the low-energy elastic properties of any hexag- 
onal [15j spontaneous spin-triplet vortex lattice, regard- 
less of the precise superconducting order parameter sym- 
metry, map onto those of the corresponding spin-singlet 
problem. The reason is that the more complicated order 
parameter, while allowing for more modes, also allows 
for additional couplings among them, which renders the 
additional modes massive. Our second main conclusion 
is that the very unusual elastic properties of the sponta- 
neous vortex lattice have easily observable consequences. 
Specifically, we predict that in ultraclean samples, the 
magnetic induction B depends nonanalytically on an ex- 



ternal magnetic field H, namely, B(H) = fiH + cH 3 / 2 , 
with [i the (linear) magnetic permeability, and c a con- 
stant. The i? 3 / 2 term is the leading nonanalyticity. In 
samples with quenched disorder strong enough to dom- 
inate the elastic properties of the flux lattice, but not 
strong enough to destroy either the superconductivity or 
the flux lattice, the nonanalyticity is the leading term as 
H -> 0, and given by B(H) oc H a with a « 0.72. This 
is the same result as in the spin-singlet case For 
very small magnetic fields either nonanalyticity is cut off 
by the lattice, which breaks the spatial rotational invari- 
ance, and the disorder-induced nonanalyticity is cut off 
at high fields by a field scale, or a corresponding length 
scale, that bounds the nonlinear elasticity regime. These 
results imply the dB/dH versus H curves shown in FigQ 
To derive these results, we start from a Landau- 
Ginzburg- Wilson (LGW) functional that allows for both 
ferromagnetic and spin-triplet superconducting order. 
The su per conducting order parameter is a matrix in spin 
space Qjji A a/3 (fe) = £^=i d M (fe) (07^02)^ . Here a, 
(3 are spin indices, k is the wave vector, and 01,2,3 are 



a 



CQ 

T3 



(a) ^S"^ 




T (b) 
II" 1 ^>-~_ 




1 1 



H cf 



H H, 



H„ 



H 



FIG. 1: Schematic representation of the predicted nonana- 
lytic behavior of B(H) for clean (a) and disordered (b) sys- 
tems. H c { is the crystal-field scale that cuts of the nonana- 
lyticity at small fields, and H n i is an upper cutoff scale deter- 
mined by the disorder. See the text for additional informa- 
tion. 
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the Pauli matrices. This relation expresses the isomor- 
phism between the complex symmetric matrix A and the 
complex 3- vector d(k). We will consider p-wave symme- 
try, which implies d fJl (k) = Ylj=i &m The tensor field 
d^j (x) is the appropriate order parameter for either He 3 
or a p-wave superconductor, d^j is characterized by 18 
real numbers, which explains the rich phenomenology of 
order parameter textures observed in He 3 (lif . 

The superconducting part of the LGW functional con- 



tains all scalars that can be constructed from d^j and 
the covariant derivative Dj = dj — iqAj(x) (q = 2e/hc), 
where the Aj are the components of the fluctuating vec- 
tor potential A, and indices in spin space and real space, 
respectively, must be contracted among themselves [lij. 
Up to bi-quadratic order in covariant gradients and ten- 
sor fields |l7| one finds (with summation over repeated 
indices implied) 



dx 



c« (Djd^i) (D^i)* + 4 2) {Djdyi) (Didtf)* + cf {Djdfuj) (Dtd*)* + t d d M d* M 



-ii^fi d d* d* -\-ii^d d* d d* 4- d d d* d* +n^d d* d d* 4- iS 5 ^d d* d* d 
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fit 



(la) 



The Cd, td, and u d are the coefficients of the LGW theory. 

To describe a ferromagnetic superconductor we need, 
in addition to d^j(x), a real vector field in spin space, 
M(x), which serves as the magnetic order parameter. 
The magnetic part of the LGW functional takes the form 
of a y 4 -theory for the fluctuating magnetization, and 
terms for the ma gnet ic field energy and the coupling be- 
tween M and A 

S im = J dx [c m (VM) 2 + t m M 2 + u m (M 2 ) 2 

+ — (V x A) 2 M- (V x A)]. (lb) 

87T J 

Finally, there is a direct coupling between the magnetic 
and superconducting order parameters (in addition to the 
indirect one via the vector potential). Considerations 
analogous to those that lead to Eq. (|lafl show that up to 
quartic order in the fields there are two such terms, 



Sfn 



dx 



M„ d v 



g 2 M ll M v d V id* 



(1c) 

where g\ and gi are two additional real coupling con- 
stants and is the Levi-Civita tensor. The complete 
LGW action is the sum of the terms in Eqs. (JTJ. 

This action is much more complicated than that for 
an s-wave superconductor. Nevertheless, the number of 
soft modes, and their long-wavelength effective action, 
are identical to those of the spin-singlet case. In par- 
ticular, the spontaneous flux lattice has the same long- 
wavelength properties. This surprising result, which is 
the basis for all of the predictions in this Letter, is inde- 
pendent of the exact nature of the spin-triplet supercon- 
ducting phase. In what follows, we first give a heuristic 
argument, and then a formal symmetry argument. 

To make the heuristic argument more transparent, let 
us consider a superconducting order parameter for the 



simplest non-unitary state, the so-called /3-state [lfij . 
(Analogous arguments can be made for other states.) It 
is given by a tensor product d — tp <8> 4> of a complex vec- 
tor ip in spin space, and a real unit vector <j> in orbital 
space, and the ground-state order parameter is given by 
i/> = Ao(l,i,0), <p = (0,0,1). In terms of ij) and 4>, the 
parts S sc and iSf m _ sc of the action read 

£f c = J dx t$ \tp\ 2 + u^p |t/>| 4 + ity |t/> x i/>*| 2 + c,/, \Dtjj\ 2 
+ \M 2 {^ (V0) 2 + d 4 , (V • 0) 2 ) + d+\{4>> £>)^| 2 ] , (2a) 



S. 



fm — s 



= dx 



-igiM ■ (i/> x V*) + 52 \M ■ ip\ 2 . (2b) 



The coefficients are simply related to those in Eqs. (Q. 

Compared to a spin-singlet s-wave superconductor, 
with a complex scalar order parameter field, we have 
more order parameter components, and also additional 
coupling terms. It turns out that the effects of these two 
complications effectively cancel each other, since the ad- 
ditional couplings generate masses for the extra degrees 
of freedom, leaving the number of massless degrees of 
freedom, and their long-wavelength Hamiltonian, identi- 
cal to that for the spontaneous flux lattice in the spin- 
singlet case. 

This can be seen as follows. We consider a region in pa- 
rameter space where the ground state of the LGW func- 
tional, Eqs. ©, i|lb|l is a vortex lattice [Io|: i.e., a lattice 
of parallel lines which each represent a topological "wind- 
ing" singularity of an overall phase associated with the 
complex vector xp. As in spin-singlet superconductors, 
the gauge couplings implicit in the covariant derivatives 
in Eq. Q2a|l force B = V x A to run along these vortex 
lines. Physically, this is the Meissner effect. Fluctuations 
in the direction of B are therefore massive. The B ■ M 
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coupling in Eq. (|lb|) then fixes the direction of M to also 
be parallel to the flux lines, up to massive fluctuations, 
and the M-(ip x t/>*) term in Eq. (|2bl) fixes the direction 
and phases of if), up to one free (Goldstone) phase whose 
vortex singularities are the flux lines. Finally, the direc- 
tion of <p is fixed by the gauge couplings in Eq. (|2a|) . 
which force <fi to also run along the flux lines. Hence, 
the entire order parameter structure is determined (up 
to massive fluctuations) by the positions of the vortex 
lines. The only Goldstone modes in the system are thus 
the two-component positional fluctuations u of the flux 
lines relative to some perfect reference lattice of parallel 
lines with an arbitrary orientation. 

A more general conclusion is reached by formal symme- 
try arguments. The superconducting part of the action, 
Eq. (Taj, is invariant under a group SO(3) xSO(3) xU(l) 
of rotations in spin sp ace, rotations in real space, and 
gauge transformations 16]. The couplings in Eq. ljlc|) 
are invariant under SO (3) co-rotations of the spin part 
of the tensor d^j and the vector M , which expresses the 
fact that both are elements of the same spin space. The 
magnetic part of the action, Eq. i|lb[l , is invariant only 
under co-rotations in spin space and in real space, so 
the entire action is invariant under a group SO(3)xU(l). 
In a ferromagnetic superconducting state this symmetry 
is spontaneously broken to SO (2), and the number of 
Goldstone modes is given by dim(SO(3)xU(l)/SO(2)) 
= 3 [19|. These are, two spin- wave-like modes due to 
the broken spin rotation symmetry, and one Anderson- 
Bogoliubov mode due to the broken gauge symmetry. 
The latter is rendered massive by means of the Higgs 
mechanism, and the former can be eliminated in favor 
of vortex-lattice degrees of freedom as shown in Ref. 
Notice that, although the triplet superconducting action 
is invariant under a much larger symmetry group than 
its singlet analog, the symmetry properties of the full ac- 
tion, and hence the number of Goldstone modes, are the 
same as in the singlet case, in agreement with the heuris- 
tic arguments given above. This symmetry argument is 
not tied to a particular superconducting ground state. 

We have corroborated the general arguments given 
above by expanding the action to Gaussian order about 
both the /3-state, and an order parameter appropriate 
for the Ai-phase in He 3, which is another nonunitary 
state. The number and nature of the soft modes found 
is consistent with the general arguments given above. 

The final effective action for any Heisenberg ferromag- 
netic p-wave superconductor [Tg| is therefore the same as 
for a ferromagnetic s-wave spin-singlet superconductor. 
It describes an Anderson-Bogoliubov mode by means of 
an 0(2) nonlinear sigma model for a phase 6, generalized 
spin waves by means of an 0(3) nonlinear sigma model 
for a unit 3-vector tp, and a coupling between the two 
by means of the gauge field A. With LGW coefficients c 



and a one thus has 

S e e = 



dx 
1 

"87T 



7T~ (V^Oz)) 2 



(V x A{x)f - M <p{x) ■ (V x A{x)) , (3) 



with Ao and Mq the average amplitudes of the supercon- 
ducting and magnetic order parameters, respectively. 

We now look for saddle-point solutions of this action 
that take the form of vortices i- e -i where the super- 
fluid velocity v = V(9 obeys the condition 



V x v(x) 
2^ 



E 



d T ^ S(x - r„(r)) 
dr 



t(x). (4) 



Here r n (r) is a parameterized line in R 3 representing the 
n'th vortex line. By minimizing the action with respect 
to 9{x) and A(x) we can express the saddle-point action 
in terms of the vortex line degrees of freedom t coupled 

to 0H3, 



? (0) 
'off 



^ J dxdy V(x - y) t(x) ■ t(y) 

- 2?r ° J dx dy V(x - y) t(x) ■ <p(y) 

2 f f 

2- / dx (V0(x)f+ 2ttM 2 A 2 / dx (V ■ <p(x)f . (5) 



Here A = 1/ y/ Ancq 2 A 2 , is the London penetration length, 
a = a-47rA 2 , and V(x) = (l/47rA 2 |a;|) exp(— |cc|/A) is a 
screened Coulomb potential. 

The action given by Eq. © can be analyzed as ex- 
plained in Ref. |l4 . The equilibrium state is a hexago- 
nal vortex lattice 15j described by two-dimensional lat- 
tice vectors R n = (X n ,Y n ). Fluctuations of the vortex 
lattice are described by a two-dimensional displacement 
field u{R n , z) such that the vortex lines are given by 

r n (z) = (X n + u x (R n ,z),Y n + u y (R n ,z),z) , (6) 

where we use z as the parameter of the line. After in- 
tegrating out the generalized spin waves, the fluctuation 
action to second order in the strain tensor |2l| 

u ij( x ) = Q [diUj(x) + dju t (x) - d a Ui(x)d a Uj(x)] (7) 



reads 



Sfluc — - 



dx 



K{d 2 z uf + 2fiUijU i: j + \{uu) 2 , (8) 



where the elastic constants re, [i (not to be confused with 
the magnetic permeability in Fig. and A (not to be 
confused with the penetration depth) can be expressed 
in terms of the coefficients of the LGW action, Eq. 10} . 
Unlike the flux lines in a conventional Abrikosov flux lat- 
tice, which are induced by an external magnetic field, the 
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flux lines in this system are spontaneously generated, and 
therefore do not have a preferred direction. As described 
for the singlet case in Ref. Q, this rotational invariance 
leads to an additional softness in the elasticity of the 
spontaneous flux lattice, due to the absence of the usual 
"tilt energy" term proportional to (d z u) 2 . This feature 
leads to anisotropic Gaussian it-propagators, where k z 
scales as fej_, with k = (k±, k z ) the wave vector. Power 
counting shows that the Gaussian theory is stable for 
all dimensions d > 5/2. This argument neglects rota- 
tional symmetry breaking by crystal fields, which pro- 
vide a long- wavelength cutoff to the applicability of our 
theory, as illustrated in Fig. 

The leading corrections to the Gaussian action are 
terms of the structure dj_u(d z u) 2 and (d z u) 4 , respec- 
tively. They lead to a least irrelevant operator with scale 
dimension — (d — 5/2), or to a wave vector dependence of 
the elastic constant fi(k± — 0, k z ) = M(l+const. x k 2d ~ 5 ). 
By arguments analogous to those given in Ref. [lj this 
leads to a nonanalytic strain-stress relation, and finally to 
the nonanalytic dependence of B on H discussed above 
in the context of Fig. 1(a). 

Even more unusual elastic properties result from the 
presence of quenched disorder. It was shown in Ref. ^3 
that ordinary impurities lead to a random-field term in 
£fl uc that couples linearly to d z u. As a result of the 
strong random-field effects, the Gaussian theory becomes 
unstable for all dimensions d < 7/2. A renormalization- 
group analysis showed that the elastic constants, as well 
as the variance of the random field, become singular func- 
tions of the wave number, and the corresponding expo- 
nents have been calculated to first order in an expansion 
in powers of e = 7/2 — d. The resulting non-Hookian, 
or nonlinear, elastic properties of the vortex lattice ex- 
tend up to a length scale £ n i = k 2 /w, where w is the 
variance of the disorder distribution. They lead to the 
predictions shown in Fig. H^b), which provide a way to 
experimentally observe the spontaneous flux lattice. 

An important question is whether there is a parame- 
ter range where the disorder is strong enough to lead to 
observable anomalous elastic effects, but not so strong as 
to destroy the p-wave superconductivity, which is very 
disorder sensitive. A necessary condition is £ <C I -C a, 
with £ the superconducting coherence length, i the mean- 
free path, and a the flux lattice constant. Let us con- 
sider ZrZn 2 13, where £ » 290l @. From the normal- 
state residual resistivity, p — 0.62/iOcm, we estimate 
I « 600 - 1000 A > f. The system is thus sufficiently 
clean to sustain p-wave superconductivity. On the other 
hand, from the value of the ordered moment, /x s = 0.17/iB 
per formula unit 6] , which gives rise to a relatively small 
spontaneous magnetic induction, B« 0.03 T, one esti- 
mates a fts y/4>o/B « 2, 500A. Here, 0o is the flux- 
quantum. The above condition is thus fulfilled in ZrZii2. 
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